BOUNDS ON THE COEFFICIENTS OF THE CHARACTERISTIC AND 

MINIMAL POLYNOMIALS 



JEAN- GUILL AUME DUMAS 

Abstract. This note presents absolute bounds on the size of the coefficients of the 
characteristic and minimal polynomials depending on the size of the coefficients of the 
associated matrix. Moreover, we present algorithms to compute more precise input- 
dependant bounds on these coefficients. Such bounds are e.g. useful to perform deter- 
ministic Chinese remaindering of the characteristic or minimal polynomial of an integer 
matrix. 



1. Introduction 

The Frobenius normal form of a matrix is used to test two matrices for similarity. Al- 
though the Frobenius normal form contains more in formation on the matrix than the 
characteristic polynomial, most efficient algorithms to compute it are based on compu- 
tations of characteristic polynomial (see for example [HI §9.7]). Now the Smith normal 
form of an integer matrix is useful e.g. in the computation of homology groups and its 
computation can be done via the integer minimal polynomial [2]. 

In both cases, the polynomials are computed first modulo several prime numbers and 
then only reconstructed via Chinese remaindering [H Theorem 10.25]. Thus, precise 
bounds on the integer coefficients of the integer characteristic or minimal polynomials of 
an integer matrix are used to know how many primes are sufficient to perform a Chinese 
remaindering of the modularly computed polynomials. Some bounds on the minimal 
polynomial coefficients, respectively the characteristic polynomial, have been presented 
in [2j, respectively in [lj. The aim of this note is to present sharper estimates in both 
cases. 

For both polynomials we present two kind of results: absolute estimates, useful to 
compare complexity constants, and algorithms which compute more precise estimates 
based on the properties of the input matrix discovered at runtime. Of course, the goal is 
to provide such estimates at a cost negligible when compared to that of actually computing 
the polynomials. 

2. Bound on the minors for the characteristic polynomial 

2.1. Hadamard's bound on the minors. The first bound of the characteristic poly- 
nomial coefficient uses Hadamard's bound, |c?et(A)| < v 'nB 2 , see e.g. jU Theorem 16.6], 
to show that the coefficients of the characteristic polynomial could be larger, but only 
slightly: 

Lemma 2.1. Let A 6 C nxn , with n > 4, whose coefficients are bounded in absolute value 
by B > 1. The coefficients of the characteristic polynomial Ca of A are denoted by Cj, 
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j = 0..n and \ \Ca\\oo = ?riaa;{|cj|}. Then 

ft 

log^H^IU) < - (log 2 (n) +log 2 (5 2 ) + 0.21163175) 

Proof. Observe that Cj, the j-th coefficient of the characteristic polynomial, is an alternate 
sum of all the (n — j) x (n — j) diagonal minors of A, see e.g. [3J, §111.7]. It is therefore 

bounded by F(n,j) = Q) \fjn - j)B 2 J . First note, that from the symmetry of the 
binomial coefficients we only need to explore the \n/2\ first ones, since y/(n — j)B 2<y ^ > 

y/W' 'for j< [n/2\ . 

The lemma is true for j = by Hadamard's bound. 

For j = 1 and n > 2, we set f(n) = \ (log 2 (F(n, 1)) - § log 2 (n) - (n - 1) log 2 (S)). 

tvt df 2n— 2+nln(n— 1)— 2n ln(n)+2 ln(n)— ln(n— 1) rpi ,i . r ji j • r 

Now = +7 — im io\ — — ■ lhus, the numerator of the derivative of 

an n z [n— 1) m(2) 

/(n) has two roots, one below 2 and one between 6 and 7. Also, f(n) is increasing from 
2 to the second root and decreasing afterwards. With n > 4, the maximal value of f(n) 
is therefore at n — 6, for which it is | log 2 (5) — | log 2 (6) < 0. 21163175. 

For other j's, Stirling's formula has been extended for the binomial coefficient by Stanica 
in [9], and gives Vi > 2, 
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J/ v2tt \j{n-j)\jj \n-j 
Now first -in 1 , -i — 1Q , 1 . -i < tttt, since the maximal value of the latter is 

12n 12j +1 12(n— 12n 6n+l ' 

at j = f . Therefore, log 2 ( ^'^"^ ) < log 2 (^) < -1.325. 

Then ls decreasing in j for 2 < j < [n/2\ so that its maximum is 2 (n-2) • 

/ \i / \ n ~i ( n —j) 

Consider now the rest of the approximation K(n,j) = (j) ( ^zj) \A n ~~ i)-^ 2 

We have log 2 (K(n, j)) = V log 2 (£ 2 ) + f log 2 (n) + f T(n,j), where T(n, j) = log 2 (^) + 

ilog 2 (^). Well T(n,j) is maximal for j = " 1+ ^ +4e " . We end with the fact that for 

n > 4 , T (n, - f log 2 (V2^) + J log 2 (^j) is maximal over Z for n = 16 

where it is lower than 0.2052. The latter is lower than 0. 21163175. □ 



We show the effectiveness of our bound on an example matrix: 
(2.1) 



11111 
11-1-1 -l 
i-i i -i -l 
i-i-i i -l 
i -l -l -l i 



This matrix has X 5 -5X 4 +40X 2 -80X+48 for characteristic polynomial and 80 = Q \/4 4 
is greater than Hadamard's bound 55.9, and less than our bound 80.66661. 

Note that this numerical bound improves on the one used in [51 lemma 2.1] since 
0.21163175 < 2 + log 2 (e) ~ 3.4427. While yielding the same asymptotic result, their 
bound would state e.g. that the coefficients of the characteristic polynomial of the example 
are lower than 21793. 



2.2. Locating the largest coefficient. The proof of lemma 12.11 suggests that the largest 
coefficient is to be found between the 0(y/n) last ones. In next lemma we take B into 
account in order to sharpen this localization. This gives a simple search procedure com- 
puting a more accurate bound on the fly, as soon as B is known. 



Lemma 2.2. Let A e C nxn ; with n > 4, whose coefficients are bounded in absolute value 
by B > 1. TTie characteristic polynomial of A is Ca- Then 



CO 

< max -i+vi+2is2 n 



where 5 ~ 5.418236. Moreover, the cost of computing the associated bound on the size is 



This localization improves by a factor close to ~, the localization of the largest coeffi- 
cient proposed in [Tj Lemma 4.1]. 



Proof. Consider F(n,j) = (") — j)B 2 ^ ^ for j = 2..|_§J- The numerator of the 
derivative of F with respect to j is 

n\y/(n-j)B 2n ~ 3 (2H(n - j) - 2H(j) - ln(n - j) - \n(B 2 ) - l) 

where H(k) = Yli=i 7 i s the fc-th Harmonic number. We have the bounds ln(fc) + 7 + 
2k \ - < H{k) < ln(fc) +7+ 2f } 1 from Theorem 2]. This bounds proves that F(n,j) 

1 — 7 3 

has at most one extremal value for 2 < j < |_§J • Moreover, ^ (n, |) < t4t — 1 + ln(^|2) 

is thus strictly negative, as soon as n > 4. Now let us define G(j) = 2H{n — j) — 2H(j) — 
\n(n — j) — ln(i? 2 ) — 1. Using the bounds on the Harmonic numbers, we have that 

~ — r < G(j)-ln ( ?^-l)+l+ln(&) < r % 

2n-2j + ^-2 2j + l U) V J 2 J 2n-2 3 + \ 2j + ^ - 2 

Then, on the one hand, we have that - — - . 2 1 — - — r-^r is increasing for 2 < ? < ^ so 
that its minimal value is Mj(n) = — 2 1 ■ — ^ at j = 2. Finally, Mj(n) > — ^ if we let 
n go to the infinity. 

On the other hand, - — %■. — r — — — \ — - is also increasing and therefore its maximal 

2n-2j + i 2^ + ^-2 

value is M s {n) = (w _ n ^+27X3n+i) at 3 = n l 2 - Finall y> M s(n) < j$5^j, its value at 
n = 4. 

Then, the monotonicity of G and its bounds prove that the maximal value of F(n,j) 
is found for j* between the solutions % an d j s of the two equations below: 



(2.2) In =l + ln(5 



(2.3) h(^)=l + MBV 2(77 - 4) 
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This proves in turn that 



. , -l + vT+^OTn. 
j* < max{0; — } 



-| 2(77-4) 

where 5 = 2e m=^) ~ 5.418236. 



Now for the complexity, we use the following recursive scheme to compute the bound: 

log(F(n, 0)) = f log(nE 2 ) 
>g(F(n, j + 1)) = log(^) + log(^) + ^ log(n - j - 1) - ^ l og (n - j) 

□ 



For instance, if we apply this lemma to matrix l2~7Ll we see that we just have to look at 
F{n,j) for j < =l±VmsMn w x 183 



3. Eigenvalue bound on the minimal polynomial 

For the minimal polynomial the Hadamard bound may also be used, but is too pes- 
simistic an estimate, in particular when the degree is small. Indeed, one can use Mignotte's 
bound on the minimal polynomial, as a factor of the characteristic polynomial. There, 
1 1 minpoly A | |oo < 2 d | |Ca| |oo, see [SI Theorem 4]. This yields that the bit size of the largest 
coefficient of the minimal polynomial is only d bits less than that of the characteristic 
polynomial. 

Therefore, one can rather use a bound on the eigenvalues determined by considera- 
tion e.g. of Gershgrin disks and ovals of Cassini (see [10J for more details on the regions 
containing eigenvalues, and [2, Algorithm OCB] for a blackbox algorithm efficiently com- 
puting such a bound). This gives a bound on the coefficients of the minimal polynomial 
of the form (3 d where (3 is a bound on the eigenvalues and d is the degree of the minimal 
polynomial. 

We can then use the following lemma to bound the coefficients of the minimal polyno- 
mial: 

Lemma 3.1. Let A G C nxn with its spectral radius bounded by [3 > 1. Let minpoly j4 (X) = 
Yfk=Q m i xl - Then 



Vi, \mA < 



'13 d ifd<(3 
min{y/j3d d ; \p^2 d (3 d } otherwise 

This improves the bound given in [21 Proposition 3.1] by a factor of log(<i) when d » (3. 



Proof. Expanding the minimal polynomial yields \rrii\ < ( d )f3 d 1 by e.g. [6, Theorem 
IV. §4.1]. Then, if d < (3, we bound the latter by (f/? d ~*. 

Now, when d > /3, we get the fist bound in two steps: first, for i < ~, we bound the 

binomial factor by d % and thus get ( d )(3 d ~ l < d % (3^~ % (3^ < dzftz since d > (3 ; second, for 

i > |, we bound the binomial factor by d d ~ l and thus get ( d )f3 d ~ r < d d ~ % f3 d ~ % < d^(3^. 

The second bound, when d > (3 is obtained by bounding the binomial coefficients 
by the middle one, (d), and using Stanica's bound [S] on the latter. This gives that 



For matrices of constant size entries, both (3 and d are 0(n). However, when d and/or (3 
is small relative to n (especially d) this may be a striking improvement over the Hadamard 
bound since the length of latter would be of order nlog(n) rather than d\og{(3). 

This is the case e.g. for the Homology matrices in the experiments of [2]. Indeed, for 
those, AA f , the Wishart matrix of A, has very small minimal polynomial degree and has 
some other useful properties which limit (3 (e.g. the matrix AA l is diagonally dominant). 
For example, the most difficult computation of [2], is that of the 25605 x 69235 matrix 
n4c6.bl2 which has a degree 827 minimal polynomial with eigenvalues bounded by 117. 
The refinement of lemma I3TT1 yields there a gain in size on the one of [2] of roughly 5%. 
In this case, this represents saving 23 modular projections and an hour of computation. 



4. Conclusion 



We have presented in this note bounds on the coefficient of the characteristic and 
minimal polynomials of a matrix. Moreover, we give algorithms with low complexity 
computing even sharper estimates on the fly. 

The refinements given here are only constant with regards to previous results but yield 
significant practical speed-ups. 
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